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To study the relativistic thermodynamic properties of a Fermi gas in a strong magnetic field, we 
construct the relativistic thermodynamic potential by the relativistic Fermi distribution function 
taking into account that the motion of particles in a plane perpendicular to the magnetic field is 
quantized. With this general potential at hand, we investigate all the thermodynamic quantities 
as a function of densities, temperatures and the magnetic field. We obtain a novel set of adiabatic 
equations. Having the expression of the pressure and adiabatic state equations, we determine the 
sound velocity for several cases revealing a new type of sound velocity. Finally, we disclose the 
magnetic cooling in the quantized electron Fermi gas, which is based on an adiabatic magnetization 
O ■ in contrast to the known adiabatic demagnetization. 
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INTRODUCTION 



The influence of strong or superstrong magnetic field on the thermodynamic properties of medium and the propa- 
gation of proper waves is an important issue in supernovae and neutron stars, the convective zone of the sun, the early 
O ■ prestellar period of the evolution of the universe, as well as in the contemporary problems of laser-matter interaction. 
Ginsburg was the first to suggest that a collapse of star must be accompanied by the generation of superstrong mag- 
netic field [1]. Based on the astrophysical data, the surface magnetic field of a neutron star is H ~ 10 11 — 10 13 G, and 
the internal field can reach H ~ 10 15 G or even higher [2]- [4]. As was shown by Bisnovati-Kogan [5], the presence of 
rotation of stars may increase the magnetic field by an additional factor of 10 3 — 10 4 . In such strong magnetic fields, it 
Oh' is expected that the thermodynamic properties and wave dynamics in degenerate electron gas would be quite different 
governed by the quantum effects. This is true when the characteristic energy of electron on a Landau level reaches the 

nonrelativistic limit of the electron chemical potential /j, = £p = h \ e \ H i i.e. H = Hs^r, where Hs = TlS = 4.4- 10 13 G 

is the Schwinger's magnetic field, vp — pp/m e = (3ir 2 ) 1 / 3 hnl^ 3 /m e is the speed of electrons at the Fermi surface, 
i <~| h is the Planck constant divided by 2-7T, m e is the electron rest mass, c is the speed of light in vacuum, | e | is the 
Oh' magnitude of the electron charge and n e is the density of electrons. 

In the past we have developed the relativistic thermodynamics of electron-ion or electron-positron-ion plasmas for 
' arbitrary temperatures in the presence of electromagnetic radiation [6] , [7] and studied numerically in Ref . [8] . Whereas 
\ the relativistic statistical thermodynamics of dense photon gas we have formulated recently [9]. In the present paper, 
we create the relativistic thermodynamics of magnetized quantum electron gas, which is becoming of increasing current 
interest motivated by its potential application in modern technology, e.g. metallic and semiconductor nanostructures 
- such as metallic nanoparticles, metal clusters, thin metal films, spintronics, nanotubes, quantum well and quantum 
dots, nano-plasmonic devices, quantum x-ray free-electron lasers, etc. Moreover, quantum electron gas is common 
in planetary interiors, in compact astrophysical objects (e.g., the interior of white dwarf stars, magnetospheres of 
' neutron stars and magnetars, etc.). 

Properties of linear electron oscillations in a dense Fermi plasma have been studied in the remote past [10]- [14]. 
Despite extensive theoretical efforts (for recent review see [15]) since then, there are questions and issues which has 
to be clarified. Answers to some salient questions are given in Refs.[16],[17], with a new type of quantum kinetic 
equations of the Fermi particles of various species and a general set of fluid equations. This kinetic equation for 
the Fermi electron gas was used to study the propagation of small longitudinal perturbations, deriving a quantum 
dispersion equation. Later the dispersion properties of linear oscillations of quantum electron-ion [18] , [19] and electron- 
positron-ion [20] plasmas, as well as of neutral -ff 3 [21] have been studied. The effects of the quantization of the orbital 
motion of electrons and the spin of electrons on the propagation of longitudinal waves in the Fermi gas have been also 
reported recently [22]. As is well known, the strong magnetic field in the Fermion gas leads to two magnetic effects 
[23] . Namely, they are the Pauli paramagnetism due to the spin of electrons and the Landau diamagnetism due to 
the quantization of the orbital motion of electrons. In Ref. [22], a novel dispersion relation of the longitudinal wave 
propagating along a magnetic field was derived, which exhibits the strong dependence on the magnetic field in radical 
contrast to the classical case. Note that in the classical limit the propagation of longitudinal wave along the magnetic 
field is unaffected by the magnetic field. 

We also note that the diamagnetic effect has a purely quantum nature and in the classical electron gas it is absent, 
because in a magnetic field the Lorentz's force -v x H acts on a particle in the perpendicular direction to a velocity 
v, so that it cannot produce work on the particle. Hence, its energy does not depend on the magnetic field. However, 
as was shown by Landau, the situation radically changes in the quantum mechanical theory of magnetism. The 
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point is that in a constant magnetic field the electrons, under the action of it, rotate in circular orbits in a plane 
perpendicular to the field H o (0,0, H ). Therefore, the motion of the electrons can be resolved into two parts: one 
along the magnetic field, in which the longitudinal component of energy is not changed and is infinite, i.e. the energy 

2 

of particle is not quantized, Eu — and the second in a plane perpendicular to H (the transverse component), in 
which the electron rotates in circular orbits with the cyclotron frequency and is finite, i.e. the transverse component of 
energy is quantized. Because the transverse motion of electrons is similar to the motion of linear harmonic oscillator, 
the corresponding energy spectrum consists of discrete energy levels [24] , [25] 

E±_ = {l+\)hu c , 

where I is the orbital quantum number [I = 0, 1, 2, ...) and u c = j s the cyclotron frequency of the electron. 
Thus, in the nonrelativistic case the net energy of electron in a magnetic field without taking into account its spin 

is 



2 ' 2m e 

If the particle has a spin, the intrinsic magnetic moment of the particle interacts directly with the magnetic field. 
The correct expression for the energy is obtained by adding an extra term flH corresponding to the energy of the 
magnetic moment jl in the field H. Hence, the electron energy levels a are determined in the nonrelativistic limit 
by the expression 

e l f = ^ + {2l+l + a)f3 B H , (1) 

where a is the operator the z component of which describes the spin orientation s = \5 (a — ±1), and /3g — ^m~c ^ s 
the Bohr magneton. 

Whereas the relativistic expression of the energy reads 



e l f = m e c\l + {21+ 1 + a)^% + 



m e c 2 m 2 c 2 



(2) 



From the expression (2) one sees that the energy spectrum of electrons consist of the lowest Landau level, t = 0, 
a = — 1, and pairs of degenerate levels with opposite polarization, u = 1. Thus each value with I ^ occurs twice, 
and that with I = once. Therefore, in the nonrelativistic limit e l e ' a can be rewritten as 



s^ = si = ^ + ehw c , (3) 



and the relativistic expression is 



ei^m e c\ 1 + 2-^ . (4) 



m e c 2 



The number of quantum states ^^Kj^c °^ a particle moving in a volume V and the interval dp z for any value of 



t in a new variables u = e{ = = VT+ Di + u 2 , D = 2-^%, Jl = = Jl + 2-^, is 

2V\e \ Hdp z _ VDdu 
(27rh) 2 c ~ 2^Af ' 

where A c = is the Compton wavelength. 



(5) 



II. THERMODYNAMICS OF FERMI GAS 



The study of thermodynamics of a Fermi gas of sufficiently low temperature and in a strong magnetic field is of 
fundamental significance. As the gas is compressed, the mean energy of electrons increases and when it becomes 
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comparable or larger than the rest energy of the electrons, the relativistic effects become important. In this case the 
temperature can be 10 3 -10 9 °K, and the relativistic degenerate electron gas will be formed. For instance, a supernova 
explosion happens when the inner part of a star first collapses to very large densities p w 10 12 — 10 14 g/cm 3 . In this 

case the Fermi energy (without a magnetic field) is ef — cpF — {Z^y^hcnJ* ~ !O s eV, which is much more than 
the thermal energy, even for the temperature T ~ 10 9 — 10 10 °K. 

We now discuss in details the density of electrons in the presence of a strong magnetic field. The equilibrium total 
density of particles is defined as 

1 J2 duf e (u,l) (6) 



with the Fermi distribution function 



exp|^| + l 

where = Kb \ and Kb is Boltzmann's constant. 

To evaluate the density n e from the expression (6), we shall consider a Fermi gas at the temperature limit | 
£huj c — n |» T. In this case the Fermi distribution function is in a good approximation described by the Heaviside 
step function H(fi — e l e ). 

Let us recall some purely quantum mechanical features of a macroscopic system. It is well known that there is 
an extremely high density of levels in the energy eigenvalue spectrum of a macroscopic system. We know also that 
the number of levels in a given finite range of the energy spectrum of a macroscopic system increases exponentially 
with the number of particles N in the system, and separations between levels are given by numbers of the 10 _JV . We 
can therefore conclude that in such case the spectrum is almost continuous, and a quasi-classical approximation is 
applicable. Thus, we can replace the summation in Eq.(6) by an integration (J^ max j^ ma * ^i) to obtain after a 
simple integration an expression of the density n e . The result is 

Zi ' ^-l + ^-tf-l-Df/ 2 } . (8) 



In equation (8) the first term is the contribution from the lowest Landau level (£ = 0), i.e. this term is associated 
with the Pauli paramagnetism and self-energy of particles. The second one results from the summation over all higher 
Landau levels. 

If the magnetic field is absent, then we recover the well known expression [24] 

' (M 2 -l) 3/2 , (9) 



3tt 2 A3 



whereas, if the magnetic field is very strong D > ("p 2 — 1), the sum in Eq.(6) vanishes and all electrons are at the 
main level, which means that the gas is fully polarized and spins of all particles are aligned opposite to the magnetic 
field. Hence, in this case for the density we have 



D 



2tt 2 A3 



, (10) 



which in the nonrelativistic limit yields [22] 



huj c m e p F 

n e = o — • (11) 



In the ultrarelativistic limit (/j, 3> 1) from Eq.(8) we obtain 

n e= PfO+o «. 1 2 (■ ( 12 ) 

From here it is clear that if (/I 2 — 1)<D, then we get the same in form expression as Eq.(ll) with the difference that 
now pf > m e c. If suppose that the density of electrons is constant, then from Eq.(ll) follows an important statement 
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that when the magnetic field increases the momentum along it decreases and a pancake configuration of Fermi energy 
becomes thin, and at H — > oo, pp goes to zero, and in the momentum space we obtain just flat-flatten. 

We now write the thermodynamic potential of relativistic electron gas with an arbitrary magnetic field and tem- 
perature 

= -2 7 i^l_^Tgy_ oo ^ln(l + e ^) . (13) 

Clearly the evolution of Eq.(13) is nontrivial, but it is possible to restrict ourselves to the most simple limiting cases. 
So that we consider the orders of magnitude of characteristic energy ef (the Fermi energy of the highest occupied 
state), T and Huj c of the system. 

We shall show that at some condition the thermodynamic potential contains a part which oscillates with a large 
amplitude as a function of H. In the remote past the oscillations of the magnetic susceptibility of metal single crystals 
at low temperatures were qualitatively predicted by Landau and experimentally detected by de Haas-van Alphen 
[24], [26]. For the oscillations a necessary condition is /i ^> Tiuj c > KbT. In which case all thermodynamic quantities 
must oscillate. It should be emphasized that at /i ^> huj c , the dependence of the not oscillatory part of thermodynamic 
potential on the magnetic field is very weak and the influence of the magnetic filed on all thermodynamic quantities 
can be ignored. 

In order to separate the oscillatory parts of £1 in Eq.(13), we use the Poisson formula 

°° -i />oo °° />oo 

O = V tf (*) = -tf (0) + / di (*) + 2Re V / di ^{t)e 2tM , (14) 
fco 2 Jl k=i Jo 

where 

In the expression (14) the first two terms describe the monotonous dependence of the thermodynamic potential on 
the magnetic field and one should take into account the magnetic field if Uoj c <~ /U, whereas the last terms depict the 
oscillatory parts, which become essential when /i ^> hu> c > KbT. We can, therefore, neglect the last terms in (14) at 
/i ~ huj c and calculate only the monotony parts. 

We consider the case, when the thermal energy KbT is much less than the Fermi energy Ep. In this case the 
distribution function is appreciably different from unity or zero only in a narrow range of values of the energy close 
to the limiting energy [i = ep. The width of this transition zone of the Fermi distribution is of the order of KbT. 
Thus, in this case the logarithm of the grand partition function (13) can be rewritten in the form 

"■^SjC*^) 8 "-^^-^-*)) ■ (15) 

The first integral denotes the value of f2 at absolute zero temperature and is the first term in the expansion of the 
corresponding quantity in powers of the small ratio T/Tp, where Tp — jf^ is the Fermi degeneracy temperature. 

The second term in (15) describes the temperature and magnetic field dependence of thermodynamical potential, 
including the first order (T/Tp) temperature correction. After integration over momentum we obtain 

n = n + nr = - 2 ^ g * (is) 
| \ CP HP + + \ E (vy/F-i-Di - (i + Di) m li+ ^~^ Dt ) + y © 2 E m } 



D 



where f( M ) = f| y_= ^t e z -l-Dl \ n=/l and I 

Taking into account the temperature correction, which follows from the last term of Eq.(16) in the limit \i ~ huj c >• 
KbT, we write down the part of thermodynamic potential that depends on the temperature 



n T = - 



Vm e c 2 D /KbT\ 2 i Ji 271 VP 2 -l-D 



12A3 Vm e c 2 y 
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If the external magnetic field equals zero, H = 0, then the expression (17) reduces to 

Vm r c 2 / KnT\ 2 



6A§ V m e c 2 

It should be noted that this equation is relativistic, from which follow expressions of the nonrelativistic limit [24] 

= _VKlT^ Fm T 



and the ultrarelativistic limit [24] 



where fi — cpp. 

We now examine the ultrarelativistic limit, i.e. 71 = -j^. 3> 1. Neglecting the temperature correction after 
integration over I (J2gLi ~^ Ji°° c ^)> we S e t 

Vm e c 2 D 2 f 17i 2 / £>\3/2l 



f2 m e c 2 D_ 2 f IT* 2 /, L>\ 3 / 2 



and thus 

p --T"w'' , t 1 + 3TK 1 -?n- (22) 

Obviously, this equation in the absence of the magnetic field, i.e. if = D = 0, reduces to the well known expression 
of the pressure 

P = l(37r 2 ) 1/3 7icn4/ 3 . (23) 

Expressions (21) and (22) describe cases, when 7i 2 > D or p 2 F > 2heH /c. In the opposite case, i.e. 7i 2 < £>, we 
obtain 

„ m„c 2 D 9 

P= w^- (24) 

Use of Eq.(ll) in Eq.(24) yields the relation, 

P=^^n 2 . (25) 

This formula shows that for given density n e the pressure of the electron gas decreases along with increase of the 
magnetic field as P ~ jj. 

If the frozen- in condition ^ = const is satisfied, then Eq.(25) casts as the equation of state of an ideal gas at the 
temperature T = const, i.e. P <~ n e . This indicates that the magnetic field sufficiently changes the equation of state. 

We next derive an expression for the thermodynamic potential or the pressure (P = —Sl/V) from Eq.(16) at T = 
in the nonrelativistic limit {pf <C m e c and Ji = 1 + -^s)- The result is 



m e c- 
\3/2 



(26) 



where 77 = ^ . 

Equation (26) describes the thermodynamic properties in the nonrelativistic limit with arbitrary magnetic field, 
which for H = reduces to 

P= ( 3 - 2 ) 2/3 ^ /3 (27) 
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and in the case, when rj > 1, reads 



(2m e e F fl 2 {2m e e F f' 2 hu: c 
P= 3rf £FV - 3^ ' (28) 



To lucidly express the pressure by the magnetic field, employing Eq.(ll) we rewrite Eq.(28) as 



C-eg) 2 , (29) 



where C = fj^ir- One can immediately see that there is a strong dependence between the pressure, the magnetic 
field and the density. Noting the relation E = — |f2, we can also define the total energy as a function of the magnetic 
field. 

With the thermodynamic potential ftp at hand from Eq.(17), we now define the entropy and specific heat for the 
Fermi electron gas as 

S = -(£) and Cv=T(^) . (30) 

\dT/ v> \dT)v,» v ' 

In the temperature limit K^T <C /i (for both cases: nonrelativistic, T <C s f /Kb = 4.35 • Wr 1Y n\J^ deg, and 
ultrarelativistic, T <C cp F /Ks = 9 • I0~ 1 nl' 3 deg) for the entropy we obtain 

cr VDK B T_j 1 2^-p 2 -l-D \ 

and the entropy per particle is 

_S__ DK B TJI f 1_ 2^/jl 2 - 1 -p \ 

S ~N~ 6A3m eC 2n e \ + D J ' 

In the absence of external magnetic field for the entropy we get the relativistic expression 



3\fn e m e c 2 ' [6) 

Obviously for the derivation of Eq.(31) we assumed that ~p 2 — 1 > D. Now if the inequality D > Jl 2 — 1 is valid for 
any value of the magnetic field, then for the entropy per electron we have 

^ ~P (33) 



Since for an adiabatic process S = const, from Eqs.(31), (32) and (33) we obtain a new adiabatic equations. First, 
the adiabatic equation without the magnetic field reads 



yf 1 + f3n e /3 = const , (34) 



where 8 = 



_ (37T 2 ) 2 / 3 ft 2 



From Eq.(34) follow well known expressions for nonrelativistic limit, which is T jnj 3 = const, and for ultrarela- 
dstic limit, which is T/n^J 3 ■ 

Next, a relativistic adiabatic equation for the case ~p 2 — 1 > D is derived from Eq.(31) as 



< 1 + -r 1 - 1 >= const , 

L noj c V sf ' J 



(35) 



, (2* 2 ) 2 ' 3 h 2 n 2 ' 3 

where e F - 2mAv+ 2 (1 _ v)3/2} 2/ 3 ■ 
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Finally, in the case of rather strong magnetic field, i.e., D > p 2 — 1 or Huj c > ep, we obtain a novel adiabatic 
equation 



TH< 



which yields in the nonrelativistic limit 



and in the ultrarelativistic case 



1 + Pj^ = const , (36) 



TH 2 

— — = const 



TH 

= const 



We specifically note that the temperature T, as follows from Eq.(36), is the function of volume and the magnetic 
field. The relation (36) implies that at constant density the increase of the magnetic field consequently leads to the 
temperature decrease as T <~ 1/H for ultrarelativistic case and T <~ 1/H 2 for nonrelativistic case. It should be 
emphasized that, therefore, this is the adiabatic magnetization process for cooling the Fermi electron gas to ultra- low 
temperatures. 

III. SOUND VELOCITY 

We now determine the velocity of sound in a nonrelativistic degenerate Fermi gas and show a very strong influence 
of the magnetic field on the expression of velocity. More importantly, we demonstrate that a new type of sound 
velocity arises, which we call the quantum magnetosound velocity, C s , that is 



2 d(p,s) fdP\ 
s d( P ,s) \dpJs' 

where p = m e n e . 

To calculate the velocity for various cases, we rewrite Eq.(37) in a new variables n e and T. Namely, 

c ,_ d(P,S) d(p,S) 



(37) 



d{n e ,T)' d(n e ,T) ' 

where the derivatives are written in Jacobian form. 

First, in the case of absence of the magnetic field, from Eq.(37) we get 



C 2 = V -f + ^-^. (38) 



3 6 h 



Next, we consider the case of a strong magnetic field, when the inequality Huj c > Ef is satisfied. In this case, Eq.(16) 
yields the following expression for the pressure 

f m 2 T 2 1 

P = hiu c n e I 1 + — f \ , (39) 



6?i n e pF 

where the second term is correction. 

Having Eqs. (33) and (39), for the nonrelativistic limit JI = 1 + after a simple calculation we obtain the 

expression of the quantum magnetosound velocity as 



L 2fi n e p F J 



We specifically note here that this is a novel expression for the sound velocity due to the strong magnetic field. 
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IV. MAGNETISM OF FERMI GASES 

In this section, we take into account a spin of electron s in the expression of energy (1). The corresponding magnetic 
moment is M — N ■ ft, where N is the total number of electrons and ft = ^ s is the magnetic moment of electron, 
with n=-n B = -0.927 • 10- 20 ^^. 

" gauss 

With the expressions of thermodynamic potentials at hand, we can find the mean magnetic moment per unit volume 
of the Fermi gas as 

V\dH'T,v,n 

We recall that, when hu c > ep, the thermodynamic potential has the form 

Vir 2 h 2 c 2 n 2 



Use of this expression (42) in Eq.(41) yields the relation 

ir 2 h 2 c 2 n|_ H 
H 2 H 

Hence the susceptibility is 



M= — -EH^^X-H. (43) 



2*2 2 2 

_tt ft c . . 

2 

It should be noted that the expression of susceptibility (44) indicates that in a strong magnetic field H > the 
gas becomes diamagnetic and the transition to the superconducting state occurs. To show this noting Eq.(44), we 
write the magnetic field induction B in a Fermi electron gas 

B = H + 4ttA? = (1 + 4ir X )H ■ (45) 

From Eqs.(44) and (45) one can immediately show the possibility of the existence of Meissner's effect, i.e. B = 

Wc estimate the magnetic field in which the metal goes to superconductor state (47rx + 1 = 0) for the density 
n e ~ 5 • 10 22 l/cm 3 

/ttW 2 \i/3 7 , , 

H = [ 2e n 2 J w 10 7 5 auss . (46) 

V. OSCILLATION OF THERMODYNAMIC POTENTIAL 

As was mentioned above, in 1930 Landau has theoretically predicted that the thermodynamic potential or the 
partition function contains an oscillatory terms in the variable 2^i/fiuj c = h/PbH. The same year the periodic 
behavior of susceptibility of metals at low temperatures and in a strong magnetic fields was experimentally discovered 
by de Haas- van Alphen. 

In order to investigate the oscillatory part of thermodynamic potential, we recall Eq.(13). Up to now we have 
evaluated Eq.(13) only in the case when e F <~ Tiu c ^> K B T. In this case the oscillatory part of Eq.(13) is very small 
in comparison with a smooth part. We now consider the case of the magnetic field in the range £f 3> huj c > KqT . 
These inequalities imply that the non-oscillatory part of the thermodynamic potential practically does not depend on 
the magnetic field and we, therefore, can use the above equations at H=0. 

For the oscillatory part of the thermodynamic potential we write the last part of Eq.(14) as 



°° />oo 

=£ / dm 

k=i J ° 



n = Yl^ Re y- I (U vj/,7), . (47) 



where 



r°° / 
y{t) = -hu c d Pz m(i + e— J. 

J — oo 
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We can rewrite the integral in Eq.(47) in dimensional variables in the form 

roo poo poo / 77 _ r e \ 

J M V(e)e t27Tkl = -hcj c m e c J du J di\n(l + ^-^-y %2M . 
If we now replace the variable I by e l e in these integrals, we obtain 

J dl *{ty 2M = .^p^ 2 J du j ^ de l e e l e ln(l + e^) e ^ k ^ , (48) 

where a = 

First, we integrate the second integral in Eq.(48). To this end, we employ a new variable e e e —~p= Qz, so that we 
rewrite the second integral as 

A = e- la Q / ( _ dz {JI + Qz) ln(l + e~ z ) e l2wk d . 

Since Jl 3> 6, the lower limit ^ 1+ ^ 2 ~^ — > — oo. Thus, in the exponent we can neglect Qz in comparison with JJ and 
therefore 

/oo _ 
ln(l + e- z ) e l± ^ fez . 

After integration twice by parts, we finally obtain the required result 



i r°° o z 
A = Q e-^ 2 " 1 ) -p— dz - e ll3z = Qjl e^-D 

P J-oo ( e + 1) 



n/3 



(e z + l) 2 ^ sinh7r/3 ' 

where /3 = 471 > efc . The advantage of this calculation is that the integrand in the second integral of Eq.(48) contributes 
practically only in the vicinity of the Fermi energy. 
The integration over u in Eq.(48) is simple, and gives 

du e- %au2 = . (49) 

V a 



Taking into account Eqs. (47), (48) and (49), finally we have for the oscillatory part of thermodynamic potential 

V(m e nu c f/ 2 ~ 1 cos[^(7I 2 -l)-f] 

27r 2 ^ 3 B fri k 3 / 2 sinh(2^kJlK B T/nw c ) ' l ° Uj 

This is the relativistic expression of f2 and indeed contains oscillatory terms in the variable (jl 2 — 1). 

In the nonrelativistic limit (Jl 2 — 1 = eF/m e c 2 ), from Eq.(50) follows a well known expression of i7 [24]. Whereas, 
in the ultrarelativistic case (JJ = 3> 1), for the cosine in Eq.(50) we have 

cos -£— fc - - . 

■m e h(jj c 4 

We now calculate the magnetic moment as the derivative of Eq.(50). To this end, we differentiate cosine alone to 
obtain 

M a y ml /2 K B Tm e c 2 (jI 2 -1) ~ I M^m e c 2 (jl 2 1) - f ] 

PB irh 2 (hiu c y/2 fr[ Vk sinh(27r 2 kJ!K B T/hu c ) ' 1 ' 



and hence the susceptibility is 



M_ ^ B r/3|mg /z m e c 2 (7I 2 -l) y, 1 — L ^ 4J 
Xosc H 7r?i 2 (^ c )3/ 2 ^ Vfcsinh(27r 2 fc7i^ B T/^ c ) ' 
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As one can see the susceptibility has an oscillatory contribution in the variable 

m e c 2 (ll 2 - 1) _ m e c 2 (lf - 1) 

)r > 1 or ^ '- > lo c , 

nuj c n 

which means that the susceptibility oscillates with a high frequency. We note here that for both nonrelativistic and 
relativistic cases these inequalities are satisfied. 

We also note that since we consider the case of /i ^> hu> c ~ KbT, the smooth part of the magnetic moment is 
smaller than the oscillatory one as [Tiu] c /m e c 2 (ji 2 — l)] 1 ^ 2 - 

We now derive the expression for 0, at low temperature KbT <C Uuj c , and k = 1, i.e. sinhx = x(l + ^-), where 
x = 2 ^ < bT < 1. The result is 

fl _ v m j"(^ efmKeiT cosr ^y-i) _ £ 

(2ir 2 ) 2 h 3 jLn e 1 6(huj c ) 2 1 1 hco c 4 J V ' 

Having the expression for the oscillatory part of thermodynamic potential (53), we can calculate the entropy, specific 
heat and the sound velocities. Namely, for the entropy we get 

* Hi),. - ^<^W=^ - J] ■ (54, 
Next, for the specific heat we have 

Cv=T(^)^ =S. (55) 



Finally, in the nonrelativistic limit (/i 2 — 1 = ^3 <C 1), for the adiabatic sound total velocity we obtain 

c 2 = ^_(dP\ _ "I , (^c) 3/2 gin[ 7Tp| TT 

S m e V9n e /5 3 27r(37r 2 n e ) 1 /3?i m i/2 m e ?iw c 4 
Clearly, as the susceptibility, the entropy, the specific heat and the sound velocity oscillate with a high frequency. 



VI. SUMMARY 



We have created the relativistic thermodynamics of magnetized quantum electron gas. To this end, we have 
derived the thermodynamic potential of relativistic electron gas with an arbitrary magnetic field and temperature, 
and investigated all the thermodynamic quantities as a function of density, temperature and the magnetic field. 
In particular cases, as the non-relativistic and the ultra-relativistic, expressions of the entropy, the specific heat, 
the pressure and the total energy are explicitly found. We have shown that the magnetic field sufficiently changes 
the equation of state. We have obtained a novel set of adiabatic equations. A novel adiabatic equation implies 
that at the constant density the increase of the magnetic field consequently leads to the temperature decrease as 
T <~ 1/H for ultrarelativistic case and T ~ 1/H 2 for nonrelativistic case. Therefore, we are here suggesting the 
adiabatic magnetization process for cooling the Fermi electron gas to ultra-low temperatures as an alternative to 
the known adiabatic demagnetization mechanism. We have also obtained expressions of sound velocity for several 
cases and demonstrated a very strong influence of the magnetic field. Moreover, we have revealed a new type of 
sound velocity, which we have named as the quantum magnetosound velocity. Finally, we have shown that in a 
strong magnetic field the Fermi gas becomes diamagnetic and the transition to the superconducting state takes place. 
The study of thermodynamics of a Fermi gas of sufficiently low temperature and in a strong magnetic field is of 
fundamental significance for understanding physics of compact astrophysical objects, e.g., the interior of white dwarf 
stars, magnetospheres of neutron stars and magnetars. The results of the present paper may be of substantial interest 
in connection with the applications in modern technology (e.g. metallic and semiconductor nanostructures) , new 
kinds of magnetic refrigerator designs, next generation intense laser-solid density matter experiments. 
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